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We investigate baryon-baryon (BB) interactions in the flavor SU(3) symmetric world, by using the 
3-flavor full QCD simulations. In the SU(3) limit, six independent two-baryon states exist for a 
given orbital angular momentum. We forcus on S-wave states and extract BB potentials using the 
method recently developed. We discuss the flavor dependence of the potentials, in particular, the 
strength of the repulsion at short distance. We find that the SU(3) flavor structure of the repulsive 
core is compatible with that expected from the Pauli principle for the quarks. 
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1. Introduction 

The generalized nuclear force, which includes not only the nucleon-nucleon (NN) interaction 
but also hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions, has been one of the central 
topics in hadron and nuclear physics. Experimental studies on the ordinary and hyper nuclei as well 
as the observational studies of the neutron stars and supernova explosions are deeply related to the 
physics of the generalized nuclear force 

As far as the NN interaction is concerned, significant number of scattering data have been 
accumulated and are parametrized in terms of the phase shifts or of the phenomenological potentials 
[0]. In particular, the moderate attraction at low energies and the strong repulsion at high energies 
seen in the S-wave NN scattering phase shift is one of the characteristic features of the NN force. 

On the other hand, the properties of the YN and YY interactions are hardly known even to- 
day due to the lack of high precision YN and YY scattering data. The approximate flavor SU(3) 
symmetry does not necessarily create a bridge between the NN force and the hyperon forces be- 
cause there are six independent flavor-channels for the scatterings between octet baryons. Although 
there have been previous theoretical attempts to fill the gap on the basis of the phenomenological 
quark models ^ and of the one-boson-exchange models [Q], it is much more desirable to study the 
generalized nuclear force from first principle QCD. 

This report is our exploratory attempt to unravel the similarities and differences among dif- 
ferent flavor representations of the baryon-baryon (BB) potentials on the basis of the lattice QCD 
simulations. To capture the essential part of physics, we take the finite but degenerated u,d,s quark 
masses, so that the system is in the flavor SU(3) limit. In the next section, we briefly review our 
method to study the baryon-baryon interaction in lattice QCD. In Section 3, we introduce the clas- 
sification of baryon-baryon systems in the flavor SU(3) limit. In Section 4, we show a setup of our 
numerical simulations. In Section 5, we show the baryon-baryon potentials in the SU(3) limit and 
discuss their implications. The last section is devoted to summary and outlook. 



2. Baryon-baryon potentials from lattice QCD 

The potential is a useful concept to describe elastic BB interactions at low-energies. Recently 
a method to extract the potential from lattice QCD simulations has been proposed and applied to 
the NN system ||5|, |6|, and ZN and AN systems [[7], |8[ ^ as well. In this method, the potential is 
defined from the Bethe-Salpeter (BS) amplitude 0s(r) through the Schrodinger equation as 

<t> E (7) = jd 3 7'U(r,?) ( t> E (r'), with T = — , k 2 = — - — (2.1) 

where pL and M are reduced and total masses of the system. Here U(r,r') is an energy independent 
but non-local potential. At low energies, it can be expanded by the local velocity v = V/jU as 
U(r,f) = V{r,v)S(r-r") = (V w + V NW + V NNL0 + ■ ■ -)8(r - r'), where N n LO term is of 0(v n ). 
Thus, in the leading-order (LO) approximation, we have 

2n <j) E {r) 
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Figure 1: A part of all 64 octet-baryon pairs. Eight 
pairs composed of the proton and octet baryons are 
expanded in several multiplets. 



27plet 


Hadron 


Mass [1/a] 


Mass [MeV] 


lOplet 


n,K 


0.5112(6) 


834.4(1.0) 




N, A, L, S 


1.068(2) 


1744.0(3.8) 



Table 1: Hadron masses in the present lattice QCD 
simulation (16 3 x 32 lattice and spatial volume of 
about (2.0 fm) 3 ). The lattice spacing a — 0.1209 fm 
is used to convert them in physical unit. 



It is shown in Ref . [ 10 ] that the LO potential V (r) defined above is a good approximation of U (r, r') 
at low energies between T ~ MeV and T ~ 45 MeV for the AW system in quenched QCD. The 
BS wave function <pE is extracted from the four point function, for example, 



W(t-t ,r) =Y,(0\p(t,x+r)n(t,x)n(t )p(t )\0) 



(2.3) 



which is dominated by the lowest energy state with the total energy Eq at large time separation, so 
that it is proportional to 0£ o (r) at large t — to. 



3. Baryon-baryon potentials in the flavor SU(3) limit 

In the flavor SU(3) limit, the ground states of a single baryon forms octet and decuplet. Sys- 
tems of two octet baryons are classified into six irreducible multiplets in flavor SU(3); 



10! 



(3.1) 



108027 8 10* 10 . 
symmetric anti-symmetric 

Since baryons obey the Fermi statistics, there exist six independent states for a given orbital angular 
momentum. In particular for the S-wave, we have singlet, octet and 27-plet for the spin-singlet state 
('So), while octet, anti-decuplet and decuplet for the spin-triplet state ( 3 Si). The six independent 
potentials for these S-wave states are the targets of our study. For example, we have computed 



J2j {0\BjB sBB |0) in order to extract y( 27 ) through eq.(2.2), where BB (21 ^ is a source operator of 
two octet baryons in the 27-plet representation, and BjBj is a sink operator with two baryons which 
has an overlap with the 27-plet component. 

We have to employ the source operator which belongs to the definite flavor representation in 
order to extract the corresponding BS wave functions at large time separation, since all two octet 
baryon states have the same energy in the infinite volume and possible energy differences due 
to the interactions in the finite volume are generally small. Fig |l| illustrates a part of all octet- 
baryon pairs. As shown in the figure, some of baryon pairs belong exclusively to 27-plet, anti- 
decuplet or decuplet. Numbers of wick contractions are relatively small for such operators, while 
more complicated source operators with the SU(3) Clebsh-Gordan coefficients are needed for other 
multiplets. Due to the limitation of computational resources, we consider only the potentials V^ 21 \ 
y(!0*) and for the moment, which can be extracted by these simpler source operators. 
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Figure 2: Effective mass plot of the 4-point func- 
tion for the two-baryon operator in the 27-plet. The 
horizontal lines stand for the threshold of two free 
baryons and its error band. 



Figure 3: Two baryon wave functions in 27-plet 
(red), anti-decuplet (green) decuplet (blue) represen- 
tations, normalized to one at the maximum distance. 



4. Numerical setup 

In order to calculate the potentials in 3-fiavor full QCD, we have utilized a gauge configuration 
set of Japan Lattice Data Grid(JLDG)/International Lattice Data Grid(ILDG) generated by CP- 
PACS/JLQCD Collaborations on 16 3 x 32 lattice with the RG improved Iwasaki gauge action at 
j8 = 1.83 and the non-perturbatively O(a) improved Wilson quark action at lattice spacing a = 0.129 
fm and spatial volume of about (2.0 fm) 3 [ ITT] , |l2| ]. The hopping parameters of the configuration 
set correspond to the flavor SU(3) symmetric point, K u = = K s = 0.13760. Quark propagators 
are calculated from the spatial wall source at to with the Dirichlet boundary condition in temporal 
direction at t = to ± 16 mod 32. The wall source is placed at 4 different time slices on each of 
800 gauge configurations, in order to enhance the signals, together with the average over forward 
and backward propagations in time. The average over discrete rotations of the cubic group is 
taken for the sink operator, in order to obtain A~\ state, whose lowest energy state is assumed 
to be dominated by the S-wave. The statistical data are divided into bins of the size eight and the 
jackknife prescription is adopted to estimate the statistical error. All of our numerical computations 
are carried out at KEK supercomputer system, Blue Gene/L and SRI 1000. 



5. Results 

Table [j] lists masses of the octet pseudo-scalar meson and baryon in our simulation, obtained 
by fitting corresponding 2-point functions. The mass of the octet pseudo-scalar meson is 834 MeV, 
which is much heavier than the physical % or K meson. For our present purpose, however, the 
flavor SU(3) symmetry is more important than the lighter u, d quark masses. 

Fig U shows the effective mass plot of the 4-point function for two baryon operator in the 27- 
plet representation. A plateau is observed at large time separation, t — to > 9, so that the 4-point 
function is well saturated by the two-baryon ground state. A similar plateaus are obtained for anti- 
decuplet and decuplet representation. We take t — to = 10 to extract baryon-baryon potentials in 
this report. 
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Figure 4: The 27-plet potential, V ' 27 ' (r), without en- Figure 5: Scattering phase shift computed from the 
ergy shift Tq. Schrtidinger equation with V ' 27 ', as a function of the 

laboratory energy 7] a j,. 



Fig|3| shows the BS wave function obtained from lattice QCD simulations. It is clearly seen that 
the wave functions depend on their flavor representation(27, 10 or 10*). An agreement between pn 
1 Sq state and pL + l So state, both of which are the members of the 27-plet, gives a good consistency 
check of our calculations. These wave functions suggest presences of repulsive interaction at short 
distance (around r < 0.4 fm) for all channels. The shape of these wave functions indicates that no 
bound state seems to exist for all channels, though more detailed analysis is needed in the finite 
volume for the definite conclusion. 

In Fig ||, the potential extracted from the BS wave function is plotted as a function of r in the 
case of the 27-plet (spin singlet) state. Note that here only the first term of eq.(2.2) is presented 
without the constant term 7b, which is expected to be small from Fig ||. The potential in Fig 
^corresponds to the NN potential in the SU(3) limit. We see a repulsive core at short distance 
and an attractive pocket around 0.6 fm. These qualitative features are similar to what was found in 
the NN system Assuming that the leading order term in the velocity expansion works well as 
demonstrated for the NN potentials |JT0|] , we can calculate the "physical" observables in the SU(3) 
symmetric world. For example, Fig ^| shows the scattering phase shift as a function of energy in the 
laboratory system, computed from the Schrodinger equation with the baryon mass of 1744 MeV in 
the present simulation and with the potential where Tq = — 10 MeV is included. 

The potential y( 10 *) without Tq in Fig^ corresponds to the NN 3 5i effective central potential 
in the SU(3) limit. We call it the "effective" central potential, since the mixing with the 3 D\ state 
is not separated from the true central potential in this calculation. We observe that y( 10 *) has a 
weaker repulsive core and a broader attractive pocket than which corresponds to the NN 1 So 
potential. Accordingly, it is more attractive as a whole. 

The potential W 10 ) without To in Fig ^ looks different from others: its repulsive core is higher 
and broader, and no clear attractive pocket is observed. Accordingly, it is much more repulsive 
than others. This feature can be also seen from the wave function in Fig ||. The stronger repulsion 
in the decuplet channel has been predicted from the phenomenological quark model where the 
decuplet states are semi-forbidden state due to the Fermi statistics of quarks Our lattice QCD 
result is compatible with this Pauli principle picture. This is more clearly seen in Fig 0, where the 



5 



Baryon-baryon potentials in the flavor SU(3) limit from lattice QCD 



Takashi Inoue 



V 



10*) 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
r[fm] 

Figure 6: Potential y< 10 *'(r) without T . 
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Figure 7: Potential V< 10 )(r) without T 
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Figure 8: Component BS amplitudes with Dirac in- Figure 9: Potential V^ l0 \r) without To, extracted 
dexes of baryon operator. Wave function of decuplet from data at three time slice, 
is given by sum of pL + Wqioi and pL + Wqho- 



4-point functions W a g a iw (r) with spinor indices aj3 , a'j3' are plotted. The spin triplet (flavor anti- 
symmetric) component is given by the sum Woifli + Woi,io- One can see that two pL + (decuplet) 
4-point functions tend to cancel each other, so that the decuplet wave function tends to vanish as if 
it were the totally forbidden state, while no such cancellation is observed for the anti-decuplet. The 
potential W 10 ) likely has a shallow attractive pocket around 0.8 fm, as seen more clearly in Fig || 
at time slice t — to = 8. 



6. Summary and outlook 

We have performed the 3-flavor full QCD simulation to study the general features of the 
baryon-baryon interaction in the flavor SU(3) limit. From the BS wave function measured on 
the lattice, we extracted the leading-order potentials in the S-wave, 

y(27) t y(lo») and V( 10 ), labeled 

by irreducible representations of flavor SU(3). All three potentials have a repulsive core at short 
distance in common, but its strength depends on the representation. In particular, W 10 ) is most 
repulsive among three, which may indicate the importance of the Pauli principle for quarks in 
baryon-baryon interactions at short distance. Also, both and V^ 10 ) show attractive pockets at 
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intermediate distance. We are currently working on the the potentials in the flavor singlet and octet 
channels to obtain a complete picture of the BB interaction in the SU(3) limit. 
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